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ABSTRACT 

We model the dark matter in galactic haloes with a self-gravitating atmosphere surrounding the 
galaxy. The galaxy serves to set the scales and the boundary conditions for the atmosphere. The 
atmosphere is treated as an isothermal Boltzmann gas, which at sufficiently large distances leads to flat 
rotation curves. Solutions to the dynamics are determined by two parameters, one of which is the ratio 
of the dark matter mass to the equilibrium temperature. From typical orbital speeds in haloes, any 
dark matter candidate that utilizes this mechanism to generate fiat rotation curves should have a mass 
to temperature ratio of ~ 400 eV per degree Kelvin. 
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The system of a galaxy plus halo is often approximated at zeroth order exclusively in terms of its 
dominant dark matter component. A simple model, which has been around for quite some time, assumes 
that all the dark matter particles are in thermal equilibrium and form a self-gravitating Boltzmann 
gas. pQ , [2] The dynamics of this system is governed by the Emden equation. [3] The standard spherically 
symmetric solution is obtained by imposing that the gravitational force vanishes at the origin r = 0. 
The resulting dark matter density profile has a 1 jr 2 behavior beyond a certain distance scale known as 
the King radius[T], which is generally associated with the galaxy core. This is just what is required for 
objects undergoing circular orbits in the medium to have an r-independent orbital speed. (Currently 
many more density profiles are available, such as the NFW profile 5], which are capable of producing 
approximately flat rotation curves at some distance scales. The profiles have been fit to observations 
in galactic haloes, and also to gravitational lensing effects. 0,13)13) The above approach assumes that 
the Boltzmann gas approximation is valid for all r, including distances that are small compared with 
the King radius; i.e. corresponding to being within the galaxy. The more involved thermodynamic 
behavior of the baryonic component to the galaxy, which interacts gravitationally with the gas, is 
ignored in this zeroth order treatment. Thus baryonic matter is to be treated as a perturbation in 
this approach. Complications can then arise in the short distance regime, as baryonic matter plays an 
important role there, and this in turn can effect the large distance behavior of the gas. Complications 
at short distances also plague dark matter simulations, where cusps appear in the density profile at 
r = 0j^5 Such singularities do not agree with observations. (In [S] it was shown that the cusps may be 
eliminated by taking into account the quantum statistics of fermions.) 

A more realistic picture should result if one takes into account the baryonic component in some 
manner from the outset. Its presence introduces a new distance scale, which we denote by Rq, as well 
as a mass scale Mq, which is the mass of the galaxy. Rq defines a boundary between the galactic 
inner region (with a low mass to luminosity ratio) and the halo outer region (with a high mass to 
luminosity ratio). We refer to the latter as the galactic atmosphere. (We shall discuss the relation of 
Rq to the King radius later.) Here we shall not be concerned with the details of the system for r < Re- 
in our model, the interior region only serves to fix the boundary conditions for the atmosphere, for 
which a simple Boltzmann gas description now appears more reasonable. This is because densities 
become low enough for the system to behave as a classical gas at sufficiently large distances from 
the origin. This gives an operational definition of the distance scale Rq, i.e., it is the distance scale 
beyond which the isothermal gas approximation is valid. The Emden equations can then be employed 
to describe the gas in the region. Here the boundary conditions should be imposed at r = Rq, rather 
than at the origin. The boundary conditions allow for a family of solutions to the Emden equations, 
which are parametrized by constants depending on Mq and Rq. The family of solutions contain the 
standard spherically symmetric solution discussed above, i.e., the one obtained by demanding that the 
gravitational force vanishes at the origin, as a special case. The resulting density profiles should then 
provide a better fit with observations than that associated with the standard spherically symmetric 
solution. As before, flat rotation curves result beyond a certain distance, which we denote by 1Z > Rq. 
Among the solutions is a well known one which gives a flat rotation curve for the entire domain. The 
solution has been referred to as the singular isothermal sphere because it is singular at r = if one 
assumes that the isothermal Boltzmann gas approximation is valid for all r. As stated above, such an 
assumption is suspect in view of the complicated dynamics within the galaxy, and we shall not presume 
that it holds. Here instead the solution is nonsingular because the domain is restricted to the galactic 
atmosphere, r > Rg- It results from particular values of the two parameters, one of which is the ratio 

' Temperature has been introduced in numerical simulations, and it was claimed to be an r— dependent function. [1] 
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of the dark particle mass to temperature, m/T. The associated fiat rotation curve for all r > Rg means 
that this solution serves as a crude fit to observations. From this fit we can obtain an estimate of the 
two parameters. Upon substituting typical orbital speeds of ~ 200 km/s found in the haloes, we find 
m/T ~ 400 eV per degree Kelvin of the atmosphere. This result does not allow for hot dark matter 
with particle mass below ~ 800 eV, since the ~ 2 °K relic background temperature serves as a lower 
bound for the equilibrium temperature of the halo. Aside from this remark, it is not our intention to 
favor any particular dark matter candidate, but rather to examine universal features resulting from this 
simple model of galactic haloes. 

In what follows we give a description of the model and some simple implications of the results. We 
begin with the basic thermodynamics of a self-gravitating isothermal gas. It results from the balancing 
of the kinetic pressure of the gas with the gravitational attraction to the galaxy and to itself. While 
thermal equilibrium may not be a good approximation within a galaxy, it most likely can be trusted 
to give a reasonable picture of the large-scale structure beyond the galaxy, and so we consider the 
restricted domain r > Rq. We shall neglect nonequilibrium processes that may occur near r = Rg 
(and also beyond some much larger scale where an atmosphere may no longer be viable). Here we shall 
assume a nonrelativistic ideal gas^] consisting of only one type of particl^J with mass m, moving in a 
slowly varying] gravitational potential cj>. The geometry of the galaxy should become less important as 
the distance from it increases. Therefore, it is reasonable to assume that spherical symmetry holds at 
lowest order, and so quantities of interest, such as <p arL d the average number density in the atmosphere 
n atrm are functions of only the radial variable. In general, the average number density has the form 
n a t m (r) = n(fi — m</>(r)), /i being the chemical potential. For a Boltzmann gas, n(/j,) — e fc s T /A 3 , where 
A = J ntfcs t ^ S t nc thermal wavelength and T is the temperature. The value of <f> at r = Rg is arbitrary. 
It is convenient to choose it to be zero so that the number density is n{fi) at r = Rq, and 

rw(r) = n atm {R G )e- m ^' kBT . (1) 

On the other hand, the value of ^ is fixed at r = Rg since it is the acceleration of gravity at r = Rg- 
Thus 

4>{R G ) =0 — = —2- , (2 

dr r=R G Rq 

Mq being the galactic mass. From the two boundary conditions we can then solve the Poisson equation 
for 4>(r), 

\Ur 2< f) - ^Gmn atm {RG)e- m ^)/k B T r > ^ (3) 

r z dr V dr J 

After a rescaling to the dimensionless variables x = and \ = q^ g ^ ne equation simplifies to 



Id/ 2 d x \ _<M , ,.s 

^T X K X di) =3Ke (4) 



along with boundary conditions 



§ We shall verify that the relativistic corrections are small. Also, we shall find the appropriate range for the dark 
particle mass where the approximation of a Boltzmann gas is valid and the effects of quantum statistics can be ignored. 

"In an atmosphere consisting of several components (including dark matter particles) in equilibrium, the heavier 
particles will sink toward the boundary region at r = Rq, while the long range behavior is governed by the lightest 
particles. 

II i.e., slowly varying with respect to the average interparticle distance 
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Here k and r are independent dimensionless parameters determining the solutions. They depend, among 
other things, on Rq and Mq. k is the ratio of the mean mass density of the atmosphere at r = Rq to the 
mean mass density of the galaxy pq — 4, S3 } while r involves m, T and the speed v x= \ = \J GMq/ Rq 
of objects undergoing circular orbits at r = Rg, 



mn(fi) 
Pa 



k R T 



(6) 



Quantities of interest are the speed v of objects undergoing circular orbits in the atmosphere as a 
function of the radial coordinate and the mass M atm of the atmosphere enclosed in a sphere of radius 
r > Rq. They are given respectively by 



\V x= iJ 



,dX 
dx 



and 



M atm (x) 

Mq 



= 3k 



dx'x l2 t 



(7) 



Eq. Q is the spherically symmetric Lane-Emden equation for infinite polytropic index, which is 
often referred to as the Emden equation. [3] Various techniques have been employed to obtain their solu- 
tions, and their application to galactic dynamics has been discussed. (See for example, [2].|10|.[TT].|12|. 
Axially symmetric solutions with rotation have also been obtained. |13j) The boundary conditions for 
this equation are standardly imposed at the origin, where the gravitational force is required to vanish 



dx 
dx 







x=0 



As mentioned above, the validity of the solutions for x < 1 is not so clear in application 
to galaxies. '"'This is not a concern here, since our domain is x > 1. Having the boundary conditions 
specified at x — 1, i.e., at the boundary of the galaxy and atmosphere, allows us to avoid complications 
in the interior region. In this model, the interior region only serves to generate the external gravitational 
field in the atmosphere. 

Solutions to the Emden equations generated from boundary conditions ^ at x = 1 form a larger 
family than those resulting from the boundary conditions at the origin stating that the gravitational 
force vanishes there. More specifically, solutions that can be consistently continued from x = 1 to the 



origin, where 



dx 



0, form a subset of the family of solutions satisfying (5 1, k and r are not 



independent for trie subset. To see this one can do another rescaling of the fielcT J and the coordinate, 
X — > X = x/ r an d x — > x — y^^x, removing the parameters from the differential equation (4), and 
putting them instead in the boundary conditions. Then 



— — (x 2 — 
x 2 dx V dx 



-x(£) 



with 



dx 



(8) 



(9) 



For generic values of k and r, the boundary conditions ([9]) cannot be generated by integrating the 

= and any x(0). [x(0) must be negative since 



Emden equation from the origin, starting from 



d.i 



x=0 



X grows with x.] Equivalently, if one were to impose the boundary condition ^¥ =0, along with 

x x—0 

some x(0) < 0, one could integrate the Emden equations to find the intersection of x(^) with the 
x— axis, along with the slope of the function at that point. One can then use ^ to numerically solve 
for specific values of k and t, defining the subset of solutions. The King radiuspQ, which is the radius 
where the dark matter density is half the density at r — 0, can be defined for these special values of n 



**X ls often set to zero, as well, at the origin, but this is a gauge condition which can be removed with a redefinition of 
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and t. It is rKi ng = y x y n ~[g) ' wnere ^atm(O) is the number density of the gas continued to the 
origin. From M it is n(jj,)e~^°\ and so 



iv»- <10) 

This expression is only valid for special values of n and r, determined from x(0) using Q. As a result of 
the latter, rxing/^G is a function only of x(0). In contrast, our approach makes no apriori restrictions 
on the values of k and r, and the behavior of the system for r < Rq is not assumed to be governed by 
the Emden equations. The King radius cannot be defined for generic values of n and t. 

Only one analytic solution to the Emden equations Q is known and it ispQ 

(11) 



X(x) = rlog ( 3 | 2 ^ 

It is called the 'singular' isothermal sphere because it is ill-defined at the origin (for any values of k 
and r). Rather than satisfying = 0, the acceleration of gravity diverges at the origin. On the 

x— 

other hand, (11) is mathematically and physically well defined for the domain of interest here, x > 1. 
It satisfies the boundary conditions (|5j) for the case of r = | and n = |, or equivalently, 



UaURG) = A^Rj T=vU (12) 

In this case, Q states that the orbital speed is independent of x, v = v x —\. It is interesting to note 
that this orbital speed is of the same order of magnitude as the rms speed u of the constituent particles 
of the gas. For a monatomic gas, u = y\v. From (7) one also gets a linear function for M atm '- 
M atm (x)/M G = x-l. 



Eq. (Ill is obviously not a solution for arbitrary values of r and k. However, (111 does describe 
the large distance behavior of generic solutions^] and satisfies the physically desired feature that the 
gravitational force tends to zero as x — > oo. (The potential is singular in this limit.) The orbital speed 
goes to the constant value, v — \f 7 2/rv x= \ (or zero, as in the case of no atmosphere k = 0) for large 
x. This is illustrated for various values of r and K [including the values r = i and K = 4 associated 



with the exact solution (111] in figures 1 and 2. M atm (x)/MG tends to a linear function of x with 



slope approaching 2r as x —¥ oo, which is illustrated in figures 3 and 4. The asymptotic region 



x > TZ/Rq > 1, where solutions tends to (11), depends on the values of k and r. The asymptotic 
region becomes the entire domain, x > 1, in the limit that (111 is exact, r — > g and k — > |. Since 



approximately flat rotation curves are observed in galactic haloes for all measured values of the radial 
coordinate beyond a certain distance scale (which we can identify with Rg), then solution (fiTl) is 



approximately valid for all x > 1. So r s» 1 and n « | roughly holds for all galaxies and (12) 
should provide a crude estimate for n(/x) and m/T. From the former, the mass density of the galactic 
atmosphere at the boundary x = 1 is approximately one-third the mass density of the galaxy. 

We can now check the validity of the nonrelativistic Boltzmann approximation which was used here. 
As we found earlier, the rms speed u of the constituent particles of the gas is of the same order of 
magnitude as v. A typical value of ~ 200 km/s can be used for the r-independent speed v = v x= \ in 
haloes, which serves as an order of magnitude estimate for It follows that relativistic corrections for 
the gas are small, although not completely negligible. (In this sense, any dark matter particle making 



tt For this we can do another change of coordinates and rewrite the differential equation in terms o f v = logx. Upon 
assuming that << ^ when y — > oo, one gets d * ~ 3ftcxp (2y — x(v) / T ): the solution of which is QXXp ■ 
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Figure 1: v/v x= i is plotted versus x for r 
ascending order of the slope at x = 1, with k = corresponding to no atmosphere. 
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Figure 2: v/v x= \ is plotted versus a; for k 



| and r 



1, | and 1. 



The values of r are listed in 



ascending order of the slope at x — 1. The limiting value of v/v x —i is V2r. 



up such a galactic atmosphere can be considered warm.) Our approximation of a classical Boltzmann 
gas is valid provided that the thermal wavelength is much smaller than the interparticle distance, or 



r, *2 2/3 



» 1 , (13) 



with the smallest value of the left- hand-side occurring at the boundary r = Rq. Upon using (12 1, 
this requires that to 8 / 3 >> 47rft 2 (pG/3) 2 / 3 /u^L 1 . For an order of magnitude estimate, we take pc ~ 
lO 9 Af /Kpc 3 along with v x= \ ~ 200 km/s, to get the result that the Boltzmann approximation is valid 
provided that dark particle mass is significantly greater than ~ 100 eV. Quantum statistics is therefore 
relevant for galactic atmospheres consisting of dark matter particles with m < 100 eV, however such 
hot dark matter scenarios are currently disfavored. (Similar conclusions appear in [5].) 



Upon substituting v x= \ ~ 200 km/s in (12) one ends up with an estimate for m/T of about 400 
eV per Kelvin. Since the ~ 2 °K background temperature provides a lower limit on the equilibrium 
temperature for any dark matter relic, the particle constituent making up the isothermal atmosphere 
cannot have a mass far below 800 eV, which is consistent with the case against hot dark matter particles. 
(A similar argument was used to estimate a lower bound on the sterile neutrino mass in [2].) The value 
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Figure 3: M atm / 'Mq is plotted versus x for r = | and k = |, | and |. The values of k are listed in 
ascending order of the slope at x = 1. 
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Figure 4: M atm /MG is plotted versus x for k — ^ and t = j, 5 and 1. The values of r are listed in 
ascending order of the slope at x = 1. The limiting value of the slope is 2r as x — > 00. 

of 800 eV is generally considered below the current lower bound for warm dark matter mass (see for 
example, [14j.[15j). and so galactic atmospheres should have an equilibrium temperature above the 
~ 2 °K background. 

The model discussed here can be applied to any dark matter candidate with mass greater than 
~ 800 eV. A promising candidate for warm dark matter is the sterile neutrino whose mass is thought to 
be in the KeV range. (For a small sample of articles, see [16] , [17] , [18] , pjj] , [20] , [2] , [22] . The possibility 
of neutral leptons being dark matter candidates was considered much earlier in [IS]-) For the case 
of galactic atmospheres consisting of 2 KeV particles, the associated equilibrium temperature should 
be ~ 5 °K. At the other end of the spectrum, an atmosphere consisting of cold dark matter particles 
would require an enormous temperature to sustain the desired density profile. For the example of a 100 
GeV particle one gets temperatures of the order 10 8-9 °K. At at such temperatures, baryonic matter 
could also play an important role in the galactic atmosphere. At a somewhat lower temperature of a 
few million degrees Kelvin, an atmosphere of 1 GeV protons could generate the l/r 2 density profile 
required for sustaining flat rotation curves. This result is not too dissimilar from early estimates of 
the temperature of baryonic gas in and around the galactic halo made in the days before the advent of 
nonbaryonic dark matter. [2"4"]. [25] Finally, we remark about a recent observation of a huge reservoir of 



ionized gas in the circumgalactic medium around the Milky Way radius, extending over a radius of 100 
kpc.[26|W| Its mass is reported to be comparable to the total baryonic mass in the disk of the galaxy 
and its temperature is of the order of 10 6 °K, which could then possibly sustain the 1/r 2 density profile 
necessary for fiat rotation curves. 
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